Exact and Approximate Sampling from the Stationary Distribution of a Markov Chain
نویسنده
چکیده
When simulating a physical system with discrete sates, one often would like to generate a sample from the stationary distribution of a Markov chain. This report focuses on three sampling methodologies which do not rely on explicitly computing the stationary distribution. Two of these lead to algorithms which can generate an exact sample in nite time. The third yields a sample whose distribution approximates, but is arbitrary close to, the stationary distribution from which one desires a sample. The approximate and one of the exact methodologies are illustrated with examples from statistical mechanics.
منابع مشابه
Exact Simulation and Error-Controlled Sampling via Regeneration and a Bernoulli Factory
Methods using regeneration have been used to draw approximations to the stationary distribution of Markov Chain Monte Carlo processes. We introduce an algorithm that allows exact sampling of the stationary distribution through the use of a regeneration method and a Bernoulli Factory to select samples within each regeneration cycle that are shown to be from the desired density. We demonstrate th...
متن کاملEecient Exact Sampling from the Ising Model Using Swendsen-wang
In Monte Carlo Markov Chain theory, a particle on a Markov chain is run for a 'long time' until the particle is nearly in the stationary distribution. Coupling from the past is a methodology which allows samples to be taken exactly from the stationary distribution of the chain, but which requires that we have a completely coupling chain. One chain often used for the ferromagnetic Ising model is...
متن کاملEecient Use of Exact Samples
Propp and Wilson (1996,1998) described a protocol called coupling from the past (CFTP) for exact sampling from the steady-state distribution of a Markov chain Monte Carlo (MCMC) process. In it a past time is identiied from which the paths of coupled Markov chains starting at every possible state would have coalesced into a single value by the present time; this value is then a sample from the s...
متن کاملOn the hardness of sampling independent sets beyond the tree threshold
We consider local Markov chain Monte-Carlo algorithms for sampling from the weighted distribution of independent sets with activity λ, where the weight of an independent set I is λ. A recent result has established that Gibbs sampling is rapidly mixing in sampling the distribution for graphs of maximum degree d and λ < λc(d), where λc(d) is the critical activity for uniqueness of the Gibbs measu...
متن کاملEfficient use of exact samples
Propp and Wilson (1996,1998) described a protocol called coupling from the past (CFTP) for exact sampling from the steady-state distribution of a Markov chain Monte Carlo (MCMC) process. In it a past time is identified from which the paths of coupled Markov chains starting at every possible state would have coalesced into a single value by the present time; this value is then a sample from the ...
متن کامل